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Abstract. We consider the steady-state regime describing the density profile of a dark 
matter halo, if dark matter is treated as a Bose-Einstein condensate. We first solve the 
fluid equation for “canonical” cold dark matter, obtaining a class of density profiles which 
includes the Navarro-Erenk-White profile, and which diverge at the halo core. We then solve 
numerically the equation obtained when an additional “quantum pressure” term is included 
in the computation of the density profile. The solution to this latter case is finite at the halo 
core, possibly avoiding the “cuspy halo problem” present in some cold dark matter theories. 
Within the model proposed, we predict the mass of the cold dark matter particle to be of 
the order of M^-c^ ~ 10“^^ eV, which is of the same order of magnitude as that predicted 
in ultra-light scalar cold dark matter models. Einally, we derive the differential equation 
describing perturbations in the density and the pressure of the dark matter fluid. 
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1 Introduction 

Observations on the rotational curves of spiral galaxies show that the velocities of the virial- 
ized material lying farther than the extent of the luminous matter from the galactic center 
reach a constant value [1]. Various theories aim at explaining this discrepancy between 
observations and Newton’s virial theorem, including a modification of the gravitational po¬ 
tential [2, 3] or of the Poisson equation [4, 5], conformal gravity [6, 7], and the metric skew 
tensor gravity [8, 9, 10]. 

Nowadays, the most promising way to explain the observations of the galactic rotation 
curves [11, 12] consists in postulating the existence of non-luminous (dark) matter, distributed 
in a halo which extends much farther than the luminous component of a galaxy. Further, 
this dark component is supposed to be non-relativistic (Cold Dark Matter, CDM [13, 14, 15], 
see also [16, 17]), since it is usually assumed to consist of massive particles with very low 
thermal velocities. Work on colliding galaxy clusters seem to confirm the existence of dark 
matter dominating the mass content of spiral galaxies and galaxy clusters [18, 19]. An 
indirect confirmation also comes from the success of the concordance cosmological model, 
or A Cold Dark Matter (ACDM) model, in reproducing the anisotropies observed in the 
cosmic microwave background [20, 21]. Among the most promising candidates for the CDM 
component are the Weakly Interacting Massive Particle (WIMP, [22]) or a population of 
zero-momentum axions [23, 24, 25, 26, 27]. 

However, these “canonical” CDM models usually feature problems in reproducing some 
observable properties of galaxies, most remarkably the under-abundance of small scale struc¬ 
ture with respect to what predicted from simulations (the “missing satellite” problem [28, 
29, 30]), the absence of a central density cusp rather predicted in numerical simulations (the 
“cusp” problem [31, 32, 33, 34]), and too little dense subhalos compared with theoretical 
expectations (the “too big to fail” problem [35]). In more detail, observations of both nearby 
dwarf galaxies and low surface brightness galaxies show that the density profile of the CDM 
halo at the core reaches a constant value [36, 37, 38, 39, 40, 41, 42]. In contrast, various 
N-body simulations predict that the CDM density distribution steepens at the center of the 
halo [43, 44, 45, 46, 47]. 
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Among the solutions proposed to overcome these issues, it has been suggested that dark 
matter could consist of a coherent scalar field with long range correlation, whose quanta are 
very light particles [48, 49, 50, 51, 52]. In fact, on short length scales, light scalar fields do 
not behave as perfect CDM and would inhibit cosmological structure growth [53, 54, 55]. 
This solution would consist of a viable mean to suppress low mass galaxies and provide cored 
profiles in CDM-dominated galaxies [51, 52, 56]. This peculiar form of dark matter might 
form a Bose-Einstein Condensate (BEC), described by the Gross-Pitaevskii or non-linear 
Schrodinger equation [57, 58, 59, 60, 61, 62, 63, 64, 65, 66, 67, 68, 69, 70, 71]. Alternatively, 
axions can also be modeled as a coherent BEC with small spatial gradient [72, 73, 74, 75, 76, 
77]. Recent 3D simulations of the gravitational collapse of wavelike cold dark matter ('I'DM) 
with a mass of the order of 10“^^eV show the effects on structure formation due to their large 
Jeans length [78, 79, 80], which improved over previous work on the subject [57, 58, 59, 60]. 
Other alternative models embed dark matter condensation into space-time with torsion [81]. 
A general review of the models proposed is given in Ref. [82]. The cosmological evolution of 
a BEC dark matter component has also been extensively explored [83, 84, 85, 86]. 

At galactic scales, the evolution of a self-gravitating CDM system can be described as a 
fluid following the equation of continuity and the Navier-Stokes Equation (NSE). When these 
equations are derived for a Bose-Einstein fluid, an additional “Quantum Pressure” (QP) term 
appears [59]. 

In this paper, we treat cold dark matter as a pressure fluid with a dynamics described 
by the NSE, and we derive the equations for the zeroth- and hrst-order perturbations in the 
density and pressure of “canonical” and BEC cold dark matter. For this, we assume a rotating 
halo in which the proper velocity of dark matter is treated as a first order perturbation in 
the motion. 

The paper is organized as follows. After the short review of fluid dynamics in Sec. 2, 
we discuss some popular halo models fitting numerical simulations in Sec. 3. In Sec. 4 we 
show that, at the lowest order, the dark matter density in the halo follows the Lane-Emden 
equation in a rotating frame. When quantum pressure is included, the Lane-Emden equation 
modifies as discussed in Sec. 5, and results for the halo prohle predict a hnite core. In 
Appendix 6, a generic expression for density perturbations and the proper velocity of dark 
matter is derived and will be part of future work. 

2 Equation for fluid dynamics 

Newton’s equations for a parcel of density p and proper velocity v, written in a reference 
frame with the 5 axis in the direction of increasing altitude, reads 

— = —Vp — V (/)-i-V • n -i-T. (2-1) 

at p 

Here, p is the pressure acting on the parcel, (j) is the gravitational potential, and r describes 
all additional external forces in the system, like the mean gravitational field generated by all 
nearby galaxies. In addition, II is a rank-two tensor describing the dissipative phenomena 
in the fluid, 

V-n = pV2v+(c + l) V (V-v). (2.2) 

The two constants appearing in Eq. (2.2) are known in the literature respectively as the 
dynamic viscosity p and the second viscosity coefficient C [87; 88]. The total time derivative 
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of the velocity field can be explicitly written as the sum of a partial time derivative and the 
dyadic (advection) term which introduces a non-linear component in Newton’s equation, 


dv 0v , civ „ 

— =- |-(v-V)v =-hV 

dt dt ^ ’ dt V 2 


- V X 


where n = |v| and we have defined the vorticity of the velocity held as 

= V X V. 


(2.3) 


(2.4) 


We assume that the galactic halo rotates at a constant rate fi, and we switch to the rotating 
frame by setting v —)■ v -|- x r, obtaining the NSE in the rotating frame 

9v 1 f \ 

— vx^ = — Vp — V — i^xfixT — 2r2xv-|-T/V. (2-5) 


Here, f2 x x r and 2 x v are respectively the Coriolis and the centrifugal acceleration terms. 
The NSE couples to two additional equations which express hux conservation (continuity 
equation), 

f + d(V-v) = 0, (2.6) 

and the value of the gravitational potential generated by the matter density p (Poisson 
equation), 

VV = 47rG/9. (2.7) 


In the following, we look for a solution to the set of Eqs. (2.5) and (2.6) in the steady-state 
regime. 


d'v dp 


( 2 . 8 ) 


Eurthermore, since we are treating the velocity as a hrst order term in the perturbation series, 
we neglect all advection terms. Under these conditions, Eq. (2.5) reads 


- V p -|- V + -|-rixr2xr-|-2rixv = pV^v. (2-9) 


while Eq. (2.6) in the steady-state regime is rewritten as the incompressibility relation V-v = 
0 for the DM how. 


3 Density profiles from N-body simulations 

In the following, we dehne the prohles so that to match the local halo density Pq = 0.4 GeV/cm^ [89, 
90, 91, 92, 93] (see also Ref. [94]) when the radius r equates the distance of the Sun from the 
center of the halo r© ~ 8.5 kpc. 

The ACDM model is very successful in predicting various observational features such 
as the galaxy rotation curves [11, 12], which are closely related to the nature of the CDM 
and the formation history of the halo. Observations show that the shape of rotation curves 
of CDM-dominated galaxies is universal across a wide mass range [37, 38, 39]. 

Early work on the analytical shape of the density profile [95, 96, 97, 98, 99, 100], 
showed that an isothermal profile with p ~ is to be expected in a flat Eriedmann- 
Robertson-Walker Universe. In addition, the density profile would steepen when sharper 
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spectra of the initial density perturbation are assumed. These claims were first validated in 
simulations [101, 102, 103, 104, 105]. The isothermal profile [100] is derived from balancing 
the gravitational and pressure forces, assuming that pressure is linearly dependent on the 
matter density, 

Piso{r) = pQ (^- 1 ) 

where q = rjrQ, qg = Ts/vq, and piso and respectively parametrize density and radial size 
of the Galaxy. The isothermal profile has the characteristic to remain finite when r —>• 0, 
approaching the value 

Pi,o(0) =P0 (l + g7")- (3.2) 

However, the isothermal profile is too shallow when compared with results from numerical 
simulations of non-colliding dark matter. Departures from the power-law behavior were 
first reported in Refs. [106, 107, 108]. A recent review has been given in Ref. [109]. Using 
numerical simulations, Navarro, Frenk, and White [NFW, 31, 43, 44] fitted results for CDM 
halos with mass spanning over four orders of magnitude with the function 

P^Fwir) = \ ^ (3.3) 

q \qs + qj 

thus predicting a rather dense central cusp growing with llq. 

Additional numerical simulations [46, 47] motivated Moore [30, 32] to propose the fol¬ 
lowing form for the CDM profile 


PM{r) 


p& ( qs + l y^^ 

\q^ + qj 


(3.4) 


which diverges at the core as ~ Ijq^^'^. The isothermal, NFW, and Moore profiles, together 
with other parametrizations such as the BE [110] and PISO [111] profiles, can collectively be 
described by the generalized expression [ 112 ] 


PCDM(g, qsa, 13, 7 ) = Pq &cDM{q, qs, a, (3, 7 ), 


with 


&CDM{q,qs,a,P,l) = q 


g? + l 
q? + 


(/3-7)/q! 


The function QcDM{q,qs,Oi, (3,'y) is normalized so that 


0CDM(l,g's,a,/3,7) = 1. 


(3.5) 

(3.6) 


(3.7) 


In Table 3, we have summarized the values of the parameters a, (3, and 7 for these popular 
halo models, including references. 

More recent N-body simulations [113, 114, 115, 116, 117, 118, 119] favor three-parameter 
profile models like the Einasto profile [120, 121 ], a generalization over the Sersic model [ 122 ], 


PE{r) = Pq 0£;(r), 


with 


&E{r) = exp 




(3.8) 


where and 6 are constants. The halo density at the core predicted by the Einasto profile 
is finite, with 


Pe{ 0) = Peexp 



(3.9) 
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Profile 

r^(kpc) 

a [3 j 

Reference 

Isothermal 

5.0 

2.0 2.0 0.0 

Bahcall [100] 

NFW 

20.0 

1.0 3.0 1.0 

Navarro et al. [31] 

Moore 

30.0 

1.0 3.0 1.5 

Moore et al. [30, 32] 

BE 

10.2 

1.0 3.0 0.3 

Binney and Evans [110] 

PISO 

4.0 

2.0 2.0 0.0 

de Boer et al. [Ill] 


Table 1. Values of the parameters appearing in Eq. (3.6) for various density profiles considered in 
the literature. 



Log 10 r / rsun 


Figure 1. The CDM profile from various fits to numerical simulations. Lines show the isothermal 
(dot-dashed line), NFW (dotted line), Moore (solid line), and Einasto profiles (dashed line). 


In Figure 1, we compare the plots of the isothermal (dot-dashed line), NFW (dotted 
line), Moore (solid line), and Einasto profiles (dashed line), with the normalization satisfying 
0cdm(?’o) = 1- For the profiles considered, we have set equal to the value given in Table 3, 
while the parameters for the Einasto profile are Vg = 20 kpc and 5 = 0.17 

4 Equations for the halo profile neglecting quantum pressure corrections 

We assume that the DM stream velocity is u <C D L, where L is the typical galactic length 
scale. To give a numerical example, for a period of rotation = 200 My and for a length 
scale L = 50 kpc, we obtain v <C 250km/s. For this reason, we first look at a numerical 
resolution of the set of Eqs. (2.5)-(2.7) in which we neglect the stream velocity. Under these 

similar comparison has been presented in Refs. [92, 123] 
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conditions, the set of equations describing the balance between pressure and density in a 
galactic CDM halo is 


V (j) = — Vp — fixfixr, 
o 

(4.1) 

II 

p) 

(4.2) 

V^(/> = 47rG/9. 

(4.3) 


Eq. (4.2) expresses mass conservation around an infinitesimal volume. The curl of Eq. (4.1) 
yields Vp x Vp = 0, which is a structural condition between the bulk pressure and density 
which is fulfilled by the barotropic relation 

p = p{p). (4.4) 

Combining the divergence of the NS equation, first line in Eq. (4.1), and the Poisson equation, 
third line in Eq. (4.1), yields the equation 

47rGp =-V +211^ (4.5) 

which is known in the literature as the Lane-Emden equation [124, 125, 126]. For the generic 
barotropic relation expressed in Eq. (4.4), the Lane-Emden equation requires a numerical res¬ 
olution [127, 128, 129]. An analytic solution exists when the relation is polytropic p oc 
with n = 1,2,5. When 11 = 0, Eq. (4.5) has often found applications in the study of 
collision-less systems such as globular clusters and primordial galaxies [1] . The rotating Lane- 
Emden equation in cylindrical coordinates has been discussed by Stodolkiewicz [130] and Os- 
triker [131] for the case of a non-rotating isothermal cylinder, by Schneider and Schmitz [132] 
for a generic polytropic fluid, and by Christodoulou and Kazanas [133] in the context of 
planetary formation for a linear polytropic relation p (x p (see also Refs. [134, 135]). 

Here, we consider the case in which the rotation is not neglected, assuming that both 
density and pressure of DM in the galactic disk do not depend on the azimuthal coordinate 
(p. We assume a cylindrical symmetry of the density, by setting 


p = Po0(r), 


(4.6) 


where 0 = 0(r) is a function depending on the distance from the galactic center r only, 
normalized so that 0(rQ) = 1. In cylindrical coordinates, Eq. (4.5) is expressed as 


AttG p + 


1 d 
r dr 



(4.7) 


To enforce the barotropic relation in Eq. (4.4), we assume that the dark matter BEC behaves 
as a self-interacting polytropic fluid, with pressure depending on density as [92, 100, 136, 137] 

p = U^p, (4.8) 


where [/ is a constant with dimensions of a velocity, known as the isothermal sound speed [133]. 
Here, we parametrize the relation between the bulk density and pressure as 

p = U^p = PQU^eir). (4.9) 
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Defining the galactic rotational period T = 27r/D, we introduce the parameter 


UJ 


= V2 VLt = 


374 My 


(4.10) 


whose order of magnitude of the parameter is approximately one. Defining 

Ut 


T = {A-kGPq) , and qu = 


rQ 


(4.11) 


and switching to the new variable q = rjrQ, the Lane-Emden Eq. (4.7) is rewritten as 


q dq \Q dq J 


(4.12) 


In the following, we solve Eq. (4.12) with the boundary conditions that &{q) and its derivative 
match the corresponding quantities from the Einasto profile at the solar neighborhood: 


0 (1) = 0^; (1) = 1, and 


dQ{q) 


dq 


dQE{q) 


q=l 


dq 


(4.13) 


9=1 


4.1 Non-rotating halo a; = 0 

When w = 0, the solution to the differential Eq. (4.12) reads 

2 


0(g) = q 


n—2 




where g^ is related to qu by 


Qs = Qu - y Qu - or qu = 


g" + l 

Ky/2^' 


(4.14) 


(4.15) 


The expression in Eq. (4.14) ensures that 0(g) > 0 for any value of g, and satisfies the 
normalization 0(1) = 1. The halo profile in Eq. (4.14) can be rephrased in terms of the 
function dehned in Eq. (3.6) as 


0(g, k) = &cDM{q, Qs, k,2 + k,2- k). 


(4.16) 


Notice that, for k = 1, Eq. (4.14) reduces to the NEW profile given in Eq. (3.3). We have 
thus obtained the important result that the NEW profile can be obtained as a solution to 
the non-rotating Lane-Emden equation. We search for the value of the parameter k which 
is consistent with the boundary condition for the derivative of the halo prohle, as given in 
Eq. (4.13). This yields to a transcendental expression for n, 


{2 - k){ 1 + q'^) + 2 k _ ^ 
1 + 9? 9f 


(4.17) 


which has to be solved numerically. Setting for example g^ = (20kpc)/rQ, we obtain k = 0.81, 
which is the value we adopt later in Sec. 4.2. In Fig. 2, we plot the value of k obtained from 
Eq. (4.17) as a function of g*. 

Fig. 3 shows the halo density profile, in units of Pq, as a function of the distance form 
the galactic center g = r/rQ. We set k = 0.2 (blue solid line), k = 0.5 (red line), k = 1 
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Figure 2. The parameter k as a function of qs = Ts/tq, as given in Eq. (4.17). 


(green line), and k = 2 (dark green line). For any value of k, all functions drop to zero for 
large values of q. At the halo core, the functions tends to infinity for 0 < k < 2 , with the 
steepness of the halo density near the galactic center decreases for increasing values of k. In 
the limit k = 2, the halo density profile at the halo core reaches a finite value at the galactic 
center, 

p{q = 0,^ = 2) = P0 , 

similarly to the prediction of the isothermal halo example in Eq. (3.1) where /3iso(0) = p©- 
Despite avoiding the cuspy halo problem, the result for k = 2 cannot be reconciled neither 
with observations nor with galactic simulation. 

The galactic rotational curves follow from the density profile in Eq. (4.14). In fact, 
inserting the density profile into the Poisson Eq. (2.7) written in the form 


d?(j) 1 dcj) 
dq^ q dq 



(4.18) 


and integrating twice from q' 
field 


0 to q' = q, we obtain the expression for the gravitational 


= In 


Qs 


q^ + q^s 


1/k 


(4.19) 


where 


Notice that, for large 
given in Ref. [138] 


= 202 km/s. 


(4.20) 


values of q ^ qs, the potential in Eq. (4.19) has the analytical expression 


■^rot 


In —, for q > qg. 


(4.21) 


The virial velocity is given by equating with the centrifugal force and reads 


V = Vrot 


' 1 - 


q^ + qf 


(4.22) 
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Figure 3. The halo density profile Q{q) = p/P q defined in Eq. (4.16), as a function of the parameter 
q = r/r© with the values k = 0.2 (blue solid line), k = 0.5 (red), k = 1 (green), and k = 2 (dark 
green). 


where Vt/) is the gradient of the potential expressed in cylindrical coordinates. We clearly 
obtain v ~ ?;rot for q S> Qs- In Fig. 4 we plot the virial velocity from Eq. (4.22) as a function 
of q for K = 0.2 (blue solid line), n = 0.5 (red line), k = \ (green line), and k = 2 (dark green 
line). We obtain that the rotational curves rise from the value u = 0 at the origin, to reach 
a flat shape for r > tq. 

4.2 Rotating halo a; / 0 

For u! ^ 0, we solve Eq. (4.12) numerically with the boundary conditions given in Eq. (4.13). 
Eig. 5 (left panel) shows the value of ©(<?) for cu = 0 (solid red), io = 0.05 (dark green), 
u = 0.5 (blue), and w = 5 (light green). Eor comparison, we have included the results for the 
isothermal profile (dot-dashed), the NFW prohle (dotted), Moore profile (solid) and Einasto 
prohle (dashed). 

The numerical resolution in the left panel of Eig. 5 shows an unrealistic behavior for 
large radii r ^ tq, where the halo density profile exponentially drops to the constant value 
instead of dropping to zero as for the case of the Einasto prohle, the NEW prohle, and 
the non-rotating solution. This behavior is due to the fact that, for r 3 > tq, the derivative 
term in Eq. (4.12) drops to zero as l/r^, and the equation reduces to 0 = We have 
further checked that, in this limit, expanding the exact solution 0 into a perturbation series, 
the hrst order perturbation also explains the exponentially damped oscillations observed in 
the full solution in Fig. 5. To sum up, since Q = uP' for large values of r, we have obtained 
an unrealistic result because the angular velocity Q, in Eq. (4.12) is assumed constant. To 
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Figure 4. The virial velocity (in km/s), given in Eq. (4.22), as a function of the parameter q = t/tq. 
We show results for k = 0.2 (blue solid line), k = 0.5 (red), k = 1 (green), and k = 2 (dark green). 




Logio r / rsun Log,o r / 


Figure 5. (Left) The function &{q) obtained from solving Eq. (4.12) with the radius-dependent 
angular momentum given in Eq. (4.26), and with the values w = 0 and k = 0.81 (red solid line), 
w = 0.05 (dark green), uj = 0.5 (blue), and w = 5 (light green). For comparison, we have also included 
the profiles obtained from the isothermal (dot-dashed line), NFW (dotted line), Moore(dashed line), 
and Einasto profiles (dashed line). (Right) The same figure as shown in the left panel, for an angular 
momentum depending on radius as in Eq. (4.26). 


circumvent this problem, we modify the Lane-Emden equation by assuming that the angular 
velocity depends on q as 


= f{q) Qq, and u}q = V2Qqt, 


(4.23) 
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where Oq is the value of n when q = 0, and f{q) is a generic function of q satisfying ^ 


/(O) = 1 and lim /(g) = 0. 

q^+oo 


(4.24) 


The Lane-Emden Eq. (4.12) is thus modihed as 

q dq \e dq J 


fiQ) 


(4.25) 


The results shown on the right panel of Eig. 5 are then obtained for the constant function 
/(g) = 1 . Here, we use the sample function 


/(<?) 


1 

T+7’ 


(4.26) 


which satishes the conditions in Eq. (4.24), to obtain a solution to Eq. (4.25) that fits the 
numerical results at large radii, as shown in the right panel of Fig. 5. 

At the halo core, the numerical resolution of Eq. (4.12) grows indefinitely for g —)• 0, 
for cj < 1 and /(g) = 1, and for any value of uj and /(g) as in Eq. (4.26). This behavior 
is similarly to what obtained for the NEW, Moore, and the halo prohle with a; = 0 shown 
in Fig. 3. Within the framework of Eq. (4.7), the divergence of p at the halo core has 
been addressed by Christodoulou and Kazanas [133], who suggest to use a composite model 
in which the hydrostatic solution only applies at large radii, while the halo distribution is 
truncated to a constant value at small radii. Here, we suggest a different solution which 
involves the addition of a quantum pressure term, as discussed in Section 5. 


5 Equations for the halo profile including quantum pressure corrections 


In the literature, the set of Eqs. (2.5)-(2.7) describes classical fluid dynamics. An analogous 
set of equations might be derived in the context of BEG cold dark matter starting from the 
Gross-Pitaevskii equation [59, 60, 63, 64, 65, 66, 67]. The proof considers the Schroedinger 
equation for the wave function describing a GDM particle, 




~2M^ + '/’(*■) + 0rot(r) + (?i>r,(r) + 


dF{p) 

dp 


/)(r,f). 


(5.1) 


where M-^ is the mass of the dark matter particle, ^(r) the gravitational potential satisfying 
the Poisson Eq. (2.7), the potential giving the Goriolis and centrifugal forces is 


i^rot — 2 ^ ^ ^ 


Jr]^ + 2f2 • V X r. 


the potential yielding the viscosities is (pn = —pv ■ Vv, and F{p) a generic function of the 
number density, p = From this latter expression, it follows that the wave function 

can be described as 

V^(r,f) = 

where S{r,t) is the action of the particle. Once the expression for has been sub¬ 

stituted into the steady-state regime of the Schroedinger Eq. (5.1), the real component of 

similar parametrization has been used in Ref. [133] 
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the expression obtained gives the Navier-Stokes Eq. (2.9) with an extra “quantum pressure” 
(QP) term appearing on its right-hand side [60, 61, 63]. Since VS" is the momentum of the 
particle, setting v = V5/M^ gives 

^ Vp -h V ^ Y -1-^^ +ClxClxr + 2ftxv = r]V‘^v + V , (5.2) 

where the pressure is given as a function of F{p) as 

P^P^^-F(p). (5.3) 

The dependence of the QP term on the inverse of the CDM particle mass squared is then 
predicted by the Schroedinger equation [59] . 

The additional QP term does not modify the barotropic relation p = p{p), since it 
does not appear in the curl of Eq. (5.2). At the same time, the divergence of Eq. (5.2) (the 
Lane-Emden equation) contains an additional term with respect to Eq. (4.5), 

4.Gp=-v(ivp)+2«^ + ^V=(Z^). (5,4) 

In cylindrical coordinates, assuming that density depends on the radial coordinate only, and 
switching to the adimensional variable q = rjrQ, Eq. (5.4) rewrites as 


0 + 


= ^2 A A J A 

q dq \Q dq J 2q dq \ y/Q dq 


i A 

q dq 


dVQ \ 
dq j 


(5.5) 


where /3 is an adimensional quantity given by 


/? = 


4 M2' 


(5.6) 


Solving Eq. (5.6) for the mass of the CDM particle and reinserting natural units, we find 




Thc^ _ 1.1 

VP 


eV, 


(5.7) 


a result within the mass range given in Refs. [53, 54, 55, 78, 79, 80]. In facts, a CDM particle 
with a mass in the range 10“^^-^ 10^“^ eV is required to solve the cusp problem and to suppress 
the small-scale spectrum [51]. Ref. [62] independently obtain a mass range similar to the one 
given in this work, in the context of galactic CDM BEC. We have shown the value of the 
mass of the CDM particle as a function of /3 in Fig. 6. 

Eq. (5.5) is solved numerically with the boundary conditions in Eq. (4.13), plus the 
additional requirement that the second and third derivatives of the numerical resolution 
match the corresponding quantities in the Einasto prohle at the solar neighborhood as well. 
For the numerical resolution, we have used the radius dependence of the angular momentum 
given in Eq. (4.26). The inclusion of QP modifies the solution to the Lane-Emden equation, 
as we show in Fig. 7 for the values w = 0 (top left, red), uj = 0.05 (top right, dark green), 
(u = 0.5 (bottom left, blue), w = 5 (bottom right, light green). For each panel, we plot the 
numerical solutions to Eq. (5.5) for the values P = 0 (solid line), P = 0.001 (dashed line, 
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Figure 6. The mass of the CDM particle as predicted from the model proposed as a function 
of the parameter (3, see Eq. (5.7). 


mass = 3.4 x 10“^^eV), (3 = 0.01 (dot-dashed line, mass = 1.1 x 10“^^eV), and 
/3 = 0.1 (dotted line, = 3.4 x 10“^^eV). For /3 / 0, the solution to Eq. (5.5) is hnte 
at the halo core. This result agrees with the assumption made in Ref. [65] that the BEC 
cold dark matter assumption alleviates the cuspy core problem appearing when simulating 
the evolution of dark matter cores. The stability of such halo, which in principle is not 
guaranteed [69], will be the subject of a future study. 

We fit the results for the value of the halo profile at the halo core from the numerical 
resolution to Eq. 5.5, obtaining the result in Eq. (5.8). The fit shows a kink at q = 0.0012, 
and reaches a constant value for (3 —>■ -|-oo, corresponding to a decreasing mass of the DM 
particle. 


0 ( 0 ) 


1907-h 


188.55 

10 . 1 / 3 -°- 99 , 


9 < 9 , 


(5.8) 


In fact, for /3 > 1, the mass density at the halo core reaches the value 


p(0) = 1907po « 760 GeV/cm^. 


(5.9) 


We show the fit in Fig. (8). 

In Appendix 6, we derive the expressions for density and pressure perturbations, as¬ 
suming that they are of the same order of magnitude as the free-streaming velocity |v|. The 
solution and discussion to this set of equations will be the subject for further study. 

6 Summary 

Ultra-light scalar dark matter with long-range correlations has been considered as a valid 
alternative to the ordinary WIMP paradigm. In this paper, we have discussed the halo 
density profile obtained when treating CDM as a perfect fluid. Our work corroborates the 
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Figure 7. The function Q{q) solution to Eq. (5.5), with the value w = 0 (top left, red), w = 0.05 
(top right, dark green), w = 0.5 (bottom left, blue), and w = 5 (bottom right, light green). For each 
panel, we have shown the results for /3 = 0 (solid line), /3 = 0.001 (dashed line), /3 = 0.01 (dot-dashed 
line), and /3 = 0.1 (dotted line). Also shown, for comparison, are the NEW profile (solid black line) 
and the Einasto profile (dashed black line). 



Figure 8. Value of 0(g) for g = 0 for different /3, as obtained from the numerical simulation (dot) 
and from the fit given in Eq. (5.8) (solid black line). 
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suggestion led by some authors [65] that the “cuspy” halo core, predicted in “canonical” 
CDM models, is removed when long-range correlations are included. For this, we have first 
solved the fluid equation for CDM, obtaining that the density profile diverges at the core, see 
Fig. 3. When an extra “quantum pressure” term arising from these long-range correlations 
is included into the fluid equation, the numerically-computed density profile remains finite 
at the halo core, see Fig. 7. 

Remarkably, for the parameters Pq ~ 0.4 GeV/cm^ and tq ~ 8.5 kpc, the model 
predicts the mass of the CDM particle to be k. 1.1 x 10“^^ eV/y^, see Eq. (5.7), 

where /? parametrizes the correlation length. For values of /3 ~ 0(1), Eq. (5.7) gives the 
mass scale of ultra-light scalar CDM, as shown in Eig. 6. 
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Appendix 


Perturbations in the fluid equations 

We linearize the NS Eq. (2.5) in the case where the density, pressure, and gravitational 
potential are perturbed as 


P = Po + Pi, P = Po+Pi, (?i = (/>o + 0i. 


( 6 . 1 ) 


Substituting this expansion into Eqs. (2.7), (2.8), and (5.2) gives the expression for the zero-th 
order perturbation as 


V (f)o =- V po — ft X Q X r, 

Po 

V^(j)o = iTTGpo, 


( 6 . 2 ) 

(6.3) 

(6.4) 


which coincide with the set of Eqs. (4.1)-(4.3). We also obtain the set of equations for the 
first order perturbations, 

lvp.-|vp„+V* + 2f! X V = + (V„p.+V. + , (6.5) 

V-v = 0, (6.6) 

VVi = 47rG/9i, (6.7) 

( 6 . 8 ) 


where the coefficients Vq, Vi, and V 2 are obtained from perturbing the term V^(y^)/y^, 
and are given by 


Vo 

Vi 


V 2 




(M. 

Po 



rpo V 2po/ 
1 

2po 


(6.9) 

( 6 . 10 ) 

( 6 . 11 ) 


15 - 





The curl of Eq. (6.5) results in the expression 


2V X X v = r/V^^, 


( 6 . 12 ) 


where we introduced the vorticity ^ = V x v. We parametrize the velocity in terms of three 
new adimensional functions u, v, w, depending on r only, as 


V = 


ur + V 6 + w z 


(6.13) 


Combining the three components of Eq. (6.12) with the continuity equation gives 


du u 

dr r 

f dv V 
( — + - 


\ dr 


w 


= 


0 , 

0 , 

0 . 


(6.14) 

(6.15) 

(6.16) 


A common solution to the set of Eqs. (6.14)-(6.16) that avoids a divergence at infinity is 
V = vo/r, for a constant vector vq. 

We now derive the expression for the divergence of Eq. (6.5). Using the barotropic 
relation in Eq. (4.8), the relation between pressure and density perturbations is 


Pi = Pu 


(6.17) 


or, writing the series expansion of the function 0 = ©o + ©i, where ©o is the solution to 
Eq. (5.5) and ©i a small perturbation, we obtain 

/9i = /9 o©i, and pi = U^PqQi. (6.18) 


Using the continuity Eq. (6.6), the divergence of the dissipation term is (V • v) = 0. 
Once the Poisson Eq. (6.7) is taken into account, the divergence of Eq. (6.5) is a differential 
equation for ©i. 


d^Oi 2 
dq^ q 



©; + 



2 ^ 
q ©0 


+ 2 



©" 

Uoj 

00 _ 


01 = 


= ^ {Oo 01 + Oi 0'i + O 2 0'/ + O 3 0?^ + O 4 , 

where the expressions for the coefficients Oi are 


(6.19) 


O'o = 6 (e'of (2 q 0'o - 30o) - 21 ©0 (©'o)^ ©" + 20©^ ©'^ ©"+ (6.20) 

+ 4©2 (©'')2 + 2©2 {3q ©' - 2©o) O'o" - q ©[]©(,'", (6.21) 

0[ = -6©o {e'of {2q ©' - 3©o) + 2Ql ©'' (7g ©' - 5©o) - 3g ©g ©'", (6.22) 

O' = ©2 [Iq (©'o)2 - Aq ©0 ©' 0 ' - 10©o ©[,] , (6.23) 

O' =-©[] (3g©' -4©o), (6.24) 

0'^ = qQl. (6.25) 
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